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We investigate the properties of quantum systems whose classical counterpart presents intermit- 
tency. It is shown, by using recent semiclassical techniques, that the quantum spectral correlations 
of such systems are expressed in terms of the eigenvalues of an anomalous diffusion operator. For 
certain values of the parameters leading to ballistic diffusion and 1/f noise the spectral proper- 
ties of our model show similarities with those of a disordered system at the Anderson transition. In 
Hamiltonian systems, intermittency is closely related to the presence of cantori in the classical phase 
space. We suggest, based on this relation, that our findings may be relevant for the description of 
the spectral correlations of (non-KAM) Hamiltonians with a classical phase space homogeneously 
filled by cantori. Finally we discuss the extension of our results to higher dimensions and their 
relation to Anderson models with long range hopping. 
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The quantum properties of a disordered system, 
namely, a non interacting particle in a random poten- 
tial are strongly affected by both the dimensionality of 
the space and the strength of disorder. In less than 
three dimensions the wavefunctions are localized in the 
thermodynamic limit for any amount of disorder. In 
three and higher dimensions there exists a metal insu- 
lator transition (MIT) for a critical amount of disorder. 
Thus for disorder below the critical one the wavefunc- 
tions are extended through the sample, the Hamiltonian 
is accurately approximated by a random matrix with the 
appropriate symmetry, and the spectral correlations are 
given by Wigner-Dyson (WD) statistics [1]. In the op- 
posite limit, wave functions are exponentially localized 
and the spectral correlations are described by Poisson 
statistics. A similar situation occurs in 'quantum chaos'. 
The celebrated Bohigas-Giannoni-Schmit conjecture [2] 
states that the WD statistics applies to the spectral cor- 
relations of quantum systems whose classical counterpart 
is chaotic. On the other hand, it is broadly accepted [3] 
that Poisson statistics describes the spectral correlations 
of quantum systems whose classical counterpart is inte- 
grable. 

Deviations from WD statistics due to wavefunction 
localization has been intensively investigated in recent 
years. In disordered systems, they are expressed through 
the dimensionless conductance g — Ec/A where Ec = 
fi/tc is the Thouless energy, tc is the classical time to cross 
the sample diffusively and A is the quantum mean level 
spacing. In the metallic regime g ^ oo and WD statis- 
tics applies. Non perturbative corrections due to a finite 
g ^ 1 were recently evaluated by Andreev and Altshuler 
[4] in the framework of the supersymmetry method [5]. 
They managed to express the two level spectral correla- 
tion function in terms of the spectral determinant of the 
classical diffusion operator. In the context of quantum 
chaos deviations from WD statistics are expected due to 
the non universality of short periodic orbits (here we do 



not discuss other sources of deviation as dynamical local- 
ization or a mixed phase space) . In a recent development 
[6], the two level spectral correlation function encoding 
such deviations was found to be related to the spectral 
determinant of the classical Perron-Frobenius operator 
which controls the evolution of the classical phase space. 
Unfortunately explicit results are hard to obtain since 
there is no a general recipe to compute the eigenvalues 
of this operator for a generic classical Hamiltonian. 

As disorder strength further increases <? ~ 1 localiza- 
tion effects become dominant and eventually the system 
undergoes a MIT. At the MIT, the wavefunction mo- 
ments Pq present anomalous scaling with respect to the 
sample size [7] L,Pq = J d''r|V'(r)|29 cx L-^''^i'^\ where 
Dq is a set of exponents describing the transition. Wave- 
functions with such a non trivial scaling arc said to be 
multifractal (for a review see [8]). Spectral fluctuations at 
the MIT (commonly referred as 'critical statistics [9] ) are 
intermediate between WD and Poisson statistics. Typi- 
cal features include: scale invariant spectrum [10], level 
repulsion and sub-Poisson number variance [11]. Differ- 
ent generalized random matrix model have been success- 
fully employed to describe 'critical statistics' [12]. 

A natural question to ask is whether critical statistics 
is related to any kind of classical motion. We shall show 
that, for a certain range of parameters, the spectral cor- 
relations of quantum systems whose classical counterpart 
presents intermittency [13] and 1/f noise are described 
by critical statistics. This is the main result of this work. 
We also suggest that our results may be useful to describe 
the spectral correlations of non-KAM Hamiltonians with 
a classical phase space homogeneously filled by cantori. 
The organization of the paper is as follows: Classical in- 
termittency is introduced by studying the evolution of a 
simple non linear map. In the context of Hamiltonian 
systems, we also discuss its relation with classical phase 
space structures Quantum spectral correlations associ- 
ated with classical intermittency are then investigated 
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by using the scmiclassical Andreev-Altshuler formalism 
above mentioned. 

Finally we discuss the extension of our results to higher 
dimensions and its relation with Anderson models with 
long range hopping. 



I. CLASSICAL INTERMITTENCY 

The phenomenon of intermittency is characterized by 
long periods of laminar (regular) motion interrupted by 
short irregular bursts [13]. As a simple example of a 
dynamical system with such features we investigate the 
following map on the real line [13,14]: Xn+i = f{xn), 
where f{x) verifies f{x) = —f{—x) and f{x + N) = 
f{x) + N with N an integer. With the above rules the 
map needs to be defined only in a restricted interval, 
f{xn) = (1 + e)x„ + a\xn\'' -1,0 < Xn < 1/2, with z 
and a real numbers and e ^ is a small control param- 
eter utilized to set the scale of the laminar phase. The 
laminar motion has its origin at points a;„ ~ where 
Xn+i ~ a;„ ± 1 and thus the orbit is transfcrcd to the 
same position in the neighboring cell. In a continuous 
time this corresponds to ballistic motion. Eventually the 
orbit leaves the region x ^ N and the dynamics becomes 
chaotic. The duration of the chaotic phase is typically 
much shorter than the ballistic one. We mention that 
due to imivcrsality [15], intermittency appears for any 
f{x) with a Taylor expansion for x <^ 1 given by the 
above relation. In [14] it was found that the density of 
probability of staying in the laminar phase a time t (or 
a distance r) has a power law tail. 
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where b = 2^a + e/2 and /U = ^^j. Indeed the mo- 
tion is superdiffusive [14] for 2 < < 3 and ballistic for 
1 < /i < 2. Based on the above result, Zumofen and 
Klafter [19] calculated the probability P{r,t) to be at 
location r at time t in the framework of the continuous 
random walk model, 
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where -P(.s, k) is the Fourier transform of P{x, t) and h ^ b 
up to factors of order the unity. For 2 < /i < 3, the above 
equation is a solution of the following fractional Fokker 
Planck equation (for a recent review see [20]), 



dP{x,t) b dt"- 



8t 2ai.l.-T^(^'*) = ^(^)^(*)- 
where we define the fractional derivative as in [20] . 
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The phenomenon of intermittency in Hamiltonian sys- 
tems has been related [21,22] to the presence of cantori 
[23] in phase space. A trajectory is typically trapped in 
the intricate cantori structures for times 'long' as com- 
pared with those involved in the transport through the 
pure chaotic phase space. Meiss and coworkers [21] de- 
scribed the transport through cantori in phase space as 
a random walk in a Bethe lattice. Such simplified model 
predicts a power law tail for the waiting times inside a 
cantori in agreement with numerical simulations [24] . De- 
tails of the distribution as the decay exponent are sensi- 
tive to the scaling properties of the cantori which depend 
on the considered energy and are thus non universal. 
Since waiting times in phase space lead to ballistic motion 
in real space [22], a trajectory in real space is a mixture of 
long ballistic flights governed by a power-law distribution 
eventually interrupted by random walks when the trajec- 
tory escapes from the cantori region. The above simpli- 
fied model cannot in principle describe the full complex- 
ity of typical KAM Hamiltonians with a mixed phase 
space where cantori with different scaling properties co- 
exist with pure chaotic and integrable components (for a 
recent investigation in this direction see [25]). A different 
situation occurs in Hamiltonians in which, due to the non 
analyticity of the classical potential, the KAM theorem 
does not hold. In certain cases [26], as a parameter is 
switched on, the whole classical phase space undergoes 
an abrupt transition from integrable to homogeneously 
filled with cantori. The absence of additional classical 
structures permits in this case utilize the formalism above 
explained though the details of P{r, t) may depend on 
the considered energy. Others class of systems with sim- 
ilar properties is that of pseudo- integrable billiards [16], 
the phase space is also fractal and the classical motion 
presents anomalous transport properties [17,18]. 



II. QUANTUM MANIFESTATIONS OF 
CLASSICAL INTERMITTENCY. 



We now investigate quantum manifestations of classi- 
cal intermittency by using semiclassical techniques. Our 
starting point is the study of the connected two level cor- 
relation function, 



i?2(s,5) = A2(p(e-s/2)p(e + s/2))-l 



(4) 



where p(e) is the density of states at energy e, A is the 
mean level spacing, the energy s is expressed in units of 
A and the averaging is over an ensemble for disordered 
systems and over an interval of energy for single deter- 
ministic chaotic systems. The spectral properties depend 
on the ratio g = E^j IS.. In the ergodic limit g ^ oo, 
the Hamiltonian can be accurately approximated by a 
random matrix with the appropriate symmetry and WD 
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statistics applies. For instance, for broken time rever- 
sal invariance, i?2(s,0) = —^^^3^^. Perturbative correc- 
tions {g ^ 1 and s 1) to this result are evaluated [11] 
by simple perturbation theory, 

Rr\s,g) = ^ReJ2P\er,,s) 

n 

where P{ f-m s) is the propagator of the diffusion equation 
(in our case (2)), e„ are the eigenvalues of the diffusion 
equation and n runs over the integers. As expected, the 
lowest mode n ~ reproduces the asymptotic form of the 
WD statistics. Non perturbative corrections (leading to 
the oscillating terms) to this result have been worked out 
by using the supersymmetry method [4] . This approach 
is valid only in the s 3> 1 limit except for systems with 
broken time reversal invariance where it is supposed to 
hold for any s. In this case. 



j, e„ = g\n 



(where periodic bound- 



R^{s,g) = Rl'^\s,g) + Rr{s,g) 



(5) 



R^"'{s,g) = cos{2-Ks) 



where D{s,g) is the spectral determinant associated 

to the diffusion equation. For normal diffusion 
D~'^{s,g) = n„5^o(l + %) ^^'^ = gn^- In principle 
one may ask whether the above formalism is applicable 
to the case of anomalous diffusion. It turns out that for 
/Lt > 2 this can be demonstrated by mapping it onto an 
Anderson model with long range disorder [28] (see below). 
For n < 2, although the mapping is in principle possible, 
we are not aware of a rigorous proof so our results should 
be considered in this case as a conjecture. 

The spectral correlations associated with classical in- 
termittency are now studied by using the above semi- 
classical techniques. We evaluate (5) for different /x 
and then we investigate the long range spectral corre- 
lations by the analysis of the number variance. We re- 
call that the number variance ^^(L) = (L)^ — (L^) = 

L + 2j^ ds{L — s)i?2(s, g) measures the stiffness of the 
spectrum. In the metallic regime, for eigenvalues sep- 
arated less than the Thouless energy, fluctuations are 
smaU and ^^(L) log(£) for L > 1. Beyond the 
Thouless energy spectral fluctuations get stronger and 
(L) ~ where d is the dimensionality of the space. 
For disorder strong enough eigenvalues are uncorrelated 
(Poisson statistics) and ^,"^{1) = L. At the MIT, the 
number variance is asymptotically proportional to 
(X < 1) with X being related with the multifractal scaling 
of the wavefunction moments [27] . 



A. Case I: 1< n < 2 

This case corresponds with classical ballistic diffusion 
(r^) oc [14]. The propagator of the diffusion equa- 
tion is given by (2), P(e„,s) = ^J-il^^ , D~^{s,g) = 



ary condition and assumed) and, by using (5), R2{s,g) ^ 
g(M-i) I s for s ^ 1. The parameter h in (2) is related 
to g as follows. For normal diffusion E^. ~ h/L/^ and 
g bL'^~^. However, in our case, since the diffusion is 
ballistic, ~ b'^'^-^/L and g ~ fe^^-^ for 1< /x < 2. We 
remark that the scale invariance of g may be modified by 
quantum corrections as in the case of a two dimensional 
weekly disordered conductor. The number variance be- 
haves asymptotically as ^^(L) ~ g(f^-'^) LlogL with a 
subleading linear term as at the MIT. Additional work 
on the wavefunctions is needed to fully characterize the 
quantum properties in this region. 



B. Case II: /tt = 2 

In this case, D-'^{s,g) = Un^^ + 7^) 
F sinh^(Wg) ^^'^ R2{s,g) can be explicitly evaluated, 

2 sin2(7rs) 



i?2(s,/i) = 



sinh {irhs) 



(6) 



where h = 1/g <^ 1. Remarkably, this correlation func- 
tion has been put forward as a definition of critical statis- 
tics [9] . It reproduces typical features at the MIT as level 
repulsion (for s <C 1 coincides with Wigner-Dyson statis- 
tics) and sub-Poisson number variance {T,'^(L) ^ hL for 
L ^ 1). In this case /i ~ 1/6 is also scale invariant. As 
discussed below, by mapping this case onto a Anderson 
model with long range disorder one can show that the 
wavefunctions are indeed multifractal as at the MIT. Fi- 
nally, we remark that the classical motion associated to 
/i = 2 leads to 1// noise [22] . 



C. Case III: 2 < < 3 

Now the dynamics is superdiffusive but sub-ballistics, 
(r^) t'^~^ [24]. Classical anomalous diffusion is de- 
scribed by the propagator, P{em s) = jj^-^, D^^{s, g) = 

J|j^_^g(l + I2-) with e„ — g\n\^^^^ . The asymptotic be- 
havior of R2{s,g) ~ §-2+1/(^-1) is power low instead of 
exponential. The conductance g is scale dependent and 
decreases as the system size increases, Ec ^ L^^/(.^^t^)^ 
g ~ I,i-2/(4-A')_ g, the power law tail of R2{s,g) 

leads to ^^(L) ~ iVC^^-i). Both the scaling of g and the 
spectral properties resemble those of a weakly disordered 
conductor in d = [28] dimensions. Finally for /x = 3 

one recovers the expected behavior S^(i) ^ xv^ of a ID 
weakly disordered conductor (no anomalous diffusion). 

We mention that similar findings have been reported 
in Anderson models with long range hopping [28]. 
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For 2 < /i < 3 the classical transport is also described by 
(2) provided that the disorder is weak enough [29]. By 
using the supersymmetry method not only the spectral 
correlations but also the wavcfunctions can be studied an- 
alytically. It turns out that the eigenfunctions are power 
law localized [28] ^ r~^~^^. In the thermodynamic 

limit, they become localized for /i > 2, delocalizcd for 
jjL <2 and multifractal for = 2 as at the MIT. Finally, 
we point out that typical features of these long range 
hopping models as power law localization and criticality 
also appears in higher dimensions [28]. Thus, in d di- 
mensions, wavefunctions are power-law localized for any 
exponent /i and multifractal for /i = c?+ 1. 

We now discuss under what conditions the above find- 
ings are relevant for deterministic Hamiltonians. Obvi- 
ously a first constraint is that the classical dynamics be 
described by (2). As discussed previously this could be 
the case for non KAM systems [26] with a classical phase 
space homogeneously filled by cantori. Although the pa- 
rameters defining (2) may depend on the considered en- 
ergy, numerical results [22] suggest that the in certain 
cases they are barely modified in a broad window of en- 
ergy and are close to the ones leading to 1// noise. From 
a quantum mechanical point view our results are only 
valid in the limit g ^ 1 where interference (not tunnel- 
ing) is the dominant quantum feature. In the context 
of quantum chaos this scale corresponds with g ~ ^W/Ti 
where h is the Planck constant and /SW is the flux swept 
across the cantorus with less flux at a given energy in one 
iteration of the map. Our results are thus applicable for 
energies such that 9 3> 1. This limit corresponds to the 
case when quantum mechanics can resolve the classical 
barrier. For 5 ~ 1 quantum mechanics cannot resolve 
the gap and in order to cross it must tunnel through it. 
We mention that recently it has been reported [30] that 
the high energy excitations of non-KAM systems as the 
Anisotropic Kepler problem or the Coulomb billiard [26] 
are correlated according to critical statistics. It would 
be interesting to check whether the classical mechanics 
of these systems present 1// noise as predicted in this 
letter. 

In conclusion, we have investigated quantum manifes- 
tation of classical intermittency. It has been shown that 
for classical ballistic diffusion and 1// noise (a special 
case of intermittency) the spectral correlations are given 
by critical statistics as at the MIT. In other cases the 
classical motion is superdiffusive but sub-ballistic, the 
wavefunctions arc power law localized and the spectral 
correlations are similar to those of a weakly disordered 
conductor in less than two dimensions. In the context 
of Hamiltonian systems we have suggested that these re- 
sults may be relevant for the description of the spectral 
correlations of non KAM systems with classical phase 
space homogeneously filled by cantori. 

Discussions with Patricio Leboeuf, Yan Fyodorov and 
Vladimir Kravtsov are gratefully acknowledged. This 



work was supported by the European Union network 
"Mathematical aspects of quantum chaos" . 



[1] M.L. Mehta, Random Matrices (Academic Press, New 

York, 1991) 2nd edition. 
[2] O. Bohigas, M.J. Gianonni and C. Schmit, Phys. Rev. 

Lett. 52, 1 (1984) 
[3] M. V. Berry and M. Tabor, Proc. Roy. Soc. London 356, 

375 (1977). 

[4] A. V. Andrccv and B. L. Altshulor, Phys. Rev. Lett. 75, 

(1995) 902. 

[5] K.B. Efetov, Adv. Phys. 32 (1983) 53. 
[e] O.Agam, B.L. Altshuler and A.V. Andreev, Phys. Rev. 
Lett. 75, (1995) 4389; E. Bogomolny and J.P. Keating, 
Phys. Rev. Lett. 77, (1996 ) 1472. 
[7] F. Wegner, Z. Phys. B 36, 209 (1980). 
[8] M. Janssen, Phys. Rep. 295, 1 (1998). 
[9] V.E. Kravtsov and K.A. Muttalib, Phys. Rev. Lett. 79, 
1913 (1997);S. Nishigaki, Phys. Rev. E 59, (1999) 2853. 
[10] B.I. Shklovskii, B. Shapiro, B.R. Sears, P. Lambrianides 

and H.B. Shore, Phys. Rev. B 47, (1993) 11487. 
[11] B.L. Altshuler, I.K. Zharekeshev, S.A. Kotochigova and 

B.I. Shklovskii, JETP 67 (1988) 62. 
[12] M. Moshe, H. Neuberger and B. Shapiro, Phys. Rev. Lett. 
73, (1994) 1497; K.A. Muttalib, Y. Chen, M.E.H. Ismail 
and V.N. Nicopoulos, Phys. Rev. Lett. 71, (1993) 471. 
[13] P. Manneville and Y. Pomeau, Phys. Lett. 75A,1 (1979). 
[14] T. Geisel, J. Nierwetberg and A. Zacherl, Phys. Rev. 

Lett. 54, (1985) 616. 
[15] I. Procaccia and H. Schuster, Phys. Rev A 28, 1210 
(1983). 

[16] P.J. Richens and M.V. Berry, Physica D 2 495 
(1981); A. Shudo and Y. Shimizu, Phy Rev. E 47 54 
(1993). 

[17] R. Artuso, I. Guarneri and L. Rebuzzini, Chaos 10 189 
(2000). 

[18] E.B. Bogomolny, U. Gerland and C. Schmit, Phys, Rev. 

E 59 R1315 (1999). 
[19] G. Zumofen and J. Klafter, Phys. Rev. E 47, 851 (1993). 
[20] R. Metzler and J. Klafter, Phys. Rep. 339, (2000) 1. 
[21] J.D. Meiss and E. Ott, Phys. Rev. Lett. 55. (1985) 2741. 
[22] T. Gcisol, A. Zacherl and G. Radons, Phys. Rev. Lett. 

59,(1987) 2503. 
[23] R.S. Mackay, J.D. Meiss and I.C. Percival, Physica D, 13 

55 (1984). 

[24] J. Klafter and G. Zumofen, Phys. Rev. E 49, 4873 (1994). 
[25] S. Denisov, J. Klafter and M. Urbakh, Phys. Rev. E, 66, 
046217 (2002). 

[26] D. Wintgen and H. Marxer, Phy Rev. Lett. 60, 971 
(1988); B.L. Altshuler and L.S. Levitov, Phys. Rep. 288, 
487 (1997) ;B. Hu, B. Li, J. Liu and Y. Gu, Phys. Rev. 
Lett. 82 (1999) 4224. 

[27] J.T. Chalker, I.V. Lerner, and R.A. Smith, Phys. Rev. 
Lett. 77, 554 (1996). 

[28] A.D. Mirlin, Y.V. Fyodorov, ot.al. Phys. Rev. E 54, 3221 

(1996) ;L.S. Levitov, Phys. Rev. Lett. 64, 547 (1990). 
[29] H.C. Fogedby, Phys. Rev. E 50, 1657 (1994). 

[30] A.M. Garcia-Garcia and J.J.M. Verbaarshot, Phys. Rev. 
E 67, 046104 (2003). 



4 



